Phase Variation of Hadronic Amplitudes 
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The phase variation with angle of hadronic amplitudes is studied with a view to understanding the 
underlying physical quantities which control it and how well it can be determined in free space. We 
find that unitarity forces a moderately accurate determination of the phase in standard amplitude 
analyses but that the nucleon-nucleon analyses done to date do not give the phase variation needed 
to achieve a good representation of the data in multiple scattering calculations. Models are examined 
which suggest its behavior near forward angles is related to the radii of the real and absorptive parts 
of the interaction. The dependence of this phase on model parameters is such that if these radii are 
modified in the nuclear medium (in combination with the change due to the shift in energy of the 
effective amplitude in the medium) then the larger magnitudes of the phase needed to fit the data 
might be attainable, but only for negative values of the phase variation parameter. 

PACS numbers: 21. 45. Be, 21.45.-v, 25.55.Ci, 13.75. Jz 

I. INTRODUCTION 

A phenomenological form often used for the nucleon-nucleon scattering amplitude is 

ika{l - ip) _i 2 

fphiq) = —e ^""^ , a = aR + tai (1) 

where a is the total cross section and p is the ratio of the real to imaginary part of the forward amplitude. This latter 
can be measured by interference with the coulomb amplitude. The parameter qr can be extracted from the fall off of 
the cross section. At a common value at the beam momentum we will consider here (1.75 GeV/c) it is usually taken 
to be 5.6 (GeV/c)~^. The parameter a/ has often been assumed to be zero for lack of better knowledge. 

Although the variation of the phase of the nucleon-nucleon amplitude with momentum transfer had been considered 
before 0,13], in seminal papers Franco and Yin Q found that the value of aj strongly influenced multiple scattering 
in light nuclei and they were able to obtain a much better representation of the data if they treated it as a free 
parameter. The values that they found were -f 10 (GeV/c)"2 and -15 (GeV/c)"^ with either value giving a dramatic 
improvement in the agreement with the data. They considered a particle scattering from ''He, "^He, deuterium and 
"'H and found the same improvement for all targets with the same values of aj. This, and other studies of multiple 
scattering \^,5, 6, 7, 8, 9], will be discussed in Section IVll 

In an attempt to estimate reasonable theoretical values for aj , Ahmad and Alvi [l0| studied an eikonal [ll| approx- 
imation based on an effective potential and concluded that a/ should be expected to have a magnitude of the order of 
1 (GeV/c)~^ or l ess, apparently in disagreement with the previous work. 

However, Ref. [lO| considered a potential in which the real and imaginary parts had the same spatial distribution, 
a Gaussian form being taken for each component. It is generally believed that the nucleon-nucleon interaction is 
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more complicated than this and that the ranges to be associated with the different parts of the interaction (real and 
imaginary) are significantly different. 

A principal aim of this paper is to investigate the effect of assuming what is hoped are reasonable estimates for 
the forms of the potential and to understand the relationship of the parameter a/ to the geometric structure of the 
interaction. 

We will often follow the eikonal method used by Ref. (lo| . In this approximation one can write 

■ik P 

fe.{q) = ^jd^b e^'i-br(6) ^ ik J hdbUqb)T{b) (2) 

where 

Tib) = 1 - e'^^^^ and xi^) = - ^ / V{^/hm^)dz, (3) 

V{r) is a complex potential and Jq is the Bessel function of zero order. 

In order to obtain the forward angle dependence of the phase we can expand Eq. [2] in powers of q and equate the 
coefficients of and in the second order expansion of these expressions with the expansion of Eq. [T] to find 

2j,-bmdb- 

We observe that there is a symmetry that exists in these expressions. If the sign of the real part of the potential is 
changed {Vaeai — > —Vaeai) then r(6) — > T*(b) which means that p — > —p and a/ — > —aj. 

We proceed in the remainder of the paper to investigate the origin and dependencies of the phase of the strong 
scattering amplitude. In section II we introduce the relationship of the phase to the difference in interaction ranges 
with a schematic model. This model, crude as it is, gives us some indication of the interdependencies among the 
physical quantities and the underlying physics involved. 

In section III we look at the conditions imposed by unitarity on the phase and its determination from phase shift 
analyses. We investigate the accuracy to which the phase can be determined in a typical realistic (K+p) case. In 
section IV we show the results of the NN analysis by Arndt et al. . 

In section V we investigate three potential models in order to give a somewhat more realistic evaluation of the 
interdependence among the parameters and their variation with energy. Section VI treats the question of determining 
the phase variation parameter, a/ from multiple scattering. In section VII we draw conclusions from the work and 
discuss the possible changes of the phase in the medium due to nucleonic and non-nucleonic mechanisms. 

II. SCHEMATIC DELTA FUNCTION MODEL 

As a first orientation, let us consider a schematic model in which the strength of the integrand in Eq. [^is concentrated 
at points in the impact parameter variable, b. In order to see the connection between the phase and the radii of the 
real and imaginary parts of the interaction it is useful to define a slightly modified amplitude such that the imaginary 
part is equal to the total cross section 

F{q) = ^/(g) along with Gib) = AiriVib) (5) 
k 



so that Eq. [2] becomes 



Note that 



F{q) = ^ / dHe^^'^Gib). (6) 



F{0)^a{p + i). (7) 



A Unitary Limits for a Gaussian Amplitude 
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We now wish to investigate the assumption that the real and imaginary parts of the strength of G{b) are concentrated 
in different regions of impact parameters, hence we consider the simple model in which the distribution of the strength 
of G{b) is expressed by (5- functions, i.e., we take 



Gib) = GnSib - bji) + iGiSib - bi). 



Now 



F(0) = GRbji + iGibi = pa + ia- 
so that the constants Gr and G/ are determined and we can write 



G(6) = a 



Ps{b-bn)+z^5{b~bj) 
bR bi 



which, using Eq. SJ leads to 



a — aji + lai 



Pbl + tbj 
2{p + i) 



or 



1 (b^-bl 



(8) 



(9) 



(10) 



(11) 



(12) 



Thus we see that a/ is directly related to the difference in the radii of the real and imaginary parts. 

If we assume that bR is to be associated with one pion exchange, dominant in this region of low momentum transfer 
p^ . and the absorption radius with the two-pion exchange range [l4l |. or about half as big, then, with a one-pion- 
exchange range of 1.4 fm we have, 



1 



1 



-(b\ - bj^) = -(0.7^ - l.V) fm^ = -0.75 fm^ « -19.3 (GeV/c) 



(13) 



so that a/ is roughly proportional to p (for small values) with a relatively large coefficient. If we were to take the 
two radii to be equal (as was done in Ref. [i3|) the coefficient would be zero. We also see a strong correlation 
between aj and p. We will see later that these general features are present in more realistic potential models. The 
relation between the phase and impact parameters has been discussed before [l^, although at much higher energies 
(PLab >100GeV/c). 

III. UNITARY CONSTRAINTS 

We now investigate to what extent the fact that physical amplitudes have an expression in partial waves with 
coefficients which satisfy unitary constraints restricts the phase parameter aj. We first consider the model amplitude 
in Eq. [T]and then treat the general case. 

A. Unitary Limits for a Gaussian Amplitude 

In the form used in Eq. [1] there are limits on the values qr and aj can take from a unitary expansion of the 
amplitude, assuming the amplitude represents scattering of a zero-spin projectile on a zero-spin target. We can write 



ika{l - ip) 12 _ ika{l ~ ip) 



Y,i\2£+l)j,{~iae)P,{x) ^l-Y^(2l+l)[St-l)Pi{x) (14) 



47r 



-e 2 



47r 



1=0 



2ik 



where 



(15) 
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^^Rt < (20) 



and X — cos 9. Thus, we can identify 

S, = l- ^!^^i^e-'=%^J,(-^aA:2). (16) 
The fact that the absolute values of Se cannot exceed unity leads to the condition 

l>\Si\^ = l- 2^iiRe + ph) + + p'){R} + Ij) (17) 

where 
and 

Ri = i?e[/e-°'='jf (-mfc2)], h = I m\i'e- '''''' jf,{-iak^)]. (19) 
For a real (= a^), so that — 0, we can write the condition 1 171 as 

2 

where Ri is the real quantity i^e~°''^^^ ji{—iaRk'^). While this condition must hold for all partial waves, numerical 
studies indicate that the s-wave unitarity is the most likely to be violated. For this case we have the result 

l-e-^-'^^<^^, (21) 

cr(l +p^) 

which can be regarded as a constraint on or p. 

This condition is similar to, but stronger than, the constraint arising from the requirement that the integrated 
elastic cross section is less than or equal to the total cross section which is 

1 - e--'^^ < (22) 
o-(l + P ) 

With values of au satisfying the condition given by Eq. [201 we included finite values of a/ and studied numerically 
the resulting values of jS'^p. It was found that for values corresponding to large partial waves (where l^^l is nearly 
unity in any case), unitarity was violated to some (often small) extent. 

We can see in the following that, in the limit of large £, unitarity in some partial waves must be violated. Using the 
limit for large t for the spherical Bessel function 

I e JH ^a^^)-^^^^e _^ - ^ [ank e) ^ ^ y ^ (2^ + ly+i ^ (^^^ 

where 

^, = -aik^+£^. (24) 
Since Rg and are going to zero with increasing £ we can drop the last term in Eq. [T7] to get 

1 > 1 - V2^/i ^ j^^^ ^ ^^^^^ [cos x^ + p sin x^] ■ (25) 

Since for some values of £ the quantity in brackets must be negative, we see that the condition will be violated for 
some S-matrix elements. 

One important caveat is that this proof holds only for the scattering of two spin-zero particles since only in that 
case can the amplitude be written in the unitarity form we have taken. The Gaussian expression for the amplitude is 
often used for a spin-averaged amplitude which is not expressible in this form. 

Thus, for the strict respect of unitarity, the form of Eq. [1] requires that a/ = and that condition [20l holds. Of 
course, one could modify the values of Sg in any partial wave in which unitarity did not hold (i.e. so that Eq. [16] is 
no longer true) but the manner of carrying out this correction is non-unique and the functional form is (perhaps only 
slightly) different from Eq. [T] We discuss this problem from a more general point of view in the next section. 



B Unitary Limits in Amplitude Analyses 

2050 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



2040 - 



"x2030 - 



2020 




2Q j^Q I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ' ' 

-0.5-0.4-0.3-0.2-0.1-0.0 0.1 0.2 0.3 0.4 0.5 

(5aj (GeV/c)-2 

FIG. 1: Values of Xd ^ function of the deviation of the phase from the original one. The dashed line indicates the minimum 
and the dotted line is drawn at a value 4 units larger 



B. Unitary Limits in Amplitude Analyses 



A similar technique to that used in the previous section can be apphed to a more general amplitude. We are 
particularly interested here in the use of a partial-wave expansion to represent experimental data and the question of 
how well the phase is determined. 

To this end, we first consider the case, again, for spin-zero on spin-zero scattering, where one attempts to change 
the phase of the amplitude by an arbitrary function, 4>{9) to obtain an amplitude with a different phase. 



f{9) = e'*('')/(0) 



(26) 



It has often been assumed that such a phase would be undetectable in elastic scattering since it does not affect the 
measurable cross section. It does, however, modify unitarity in a manner similar to that seen in the previous section. 
We will take the form of (j)(9) to be linear in t, i.e. 



(27) 



= -i(5a/t = ifc/2fc^(l - COS0) 

where ly — 5a jk^. We can expand both amplitudes in Eq. [26] in Legendre series to find 

fiO) = ^/£P^(cos0) and /(0) = ^/lPl(cos0) 



(28) 



where 



„ Sg — I ~ Sl — ^ 

h and /l = 



2ik 



2ik 



(29) 



Using Bauer's series. 



= ^i^(2A + l)PA(cos0)jA(-z^), 

A 



(30) 
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FIG. 2; Comparison of phases of K^p scattering. Upper panel: The sohd curve gives the phase from the original amplitude 
and the dashed curve the result of modifying it with an additional phase with 5ai = —0.51 (GeV /c)~^ . The dash-dot curve 
gives the phase of the amplitude resulting to the best fit to the experimental data plus the phase data out to 0.32 (GeV/c)^. 
The phases match out to about 0.36 (GeV/c)^. The case chosen for this comparison is the worst case (has the largest value 
of from the phase "data"). Lower panel; Original phase and limits corresponding to about two standard deviations in the 
classical estimate. The values that correspond to this point on Fig. [T]are about ± 0.2 (GeV/c)~^. 



and Eqs. \Tl\ and we can express the coefficients in the expansion of f{9) as 



'0,0,0 



(31) 



where C 



Mi,M2,M3 



is a Clebsch-Gordon coefficient. Solving for the S- matrix element Sl from Eqs. [29l and [3T] one can 



check if 15*^1 < 1. We find, with numerical studies, that unitarity is always violated in some partial wave. 

How to judge the seriousness of a violation of unitarity is perhaps not obvious. One way to do so is to correct the 
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violation and see what difference the change makes in the cross section and other observables derived from the new 
amplitude. This leads then to constraints due to the data. Rather than make ad hoc changes to the S-matrix elements, 
it is preferable to perform a search fitting the data to make the decision about how the parameters are to be changed 
to preserve unitarity, give the desired phase and, at the same time, give the best fit to the data in the sense of a lowest 

. .... . n 

We have implemented this idea in a fitting program for K"*"p scattering |16| . a system for which the data are relatively 

good. The isospin structure is the same as in nuclcon-nucleon scattering but the spin structure is less complex. For 
K"'"p scattering there is no one-pion-exchange (OPE) contribution as there is in the nucleon-nucleon case. It has been 
suggested that OPE will help to determine the phase of the amplitude for NN scattering but in the present test 
we have only the unitarity constraints to determine the phase. 

We now study how well the phase can be determined in the process of finding a fit of a unitary form to the data. 
For K"'"p scattering we can write the amplitude as F(6) + a ■ nG{6) where n is a unit vector perpendicular to the 
scattering plane. In the course of this work a minimum lower than those found in Ref. was observed. The original 
best found was 2031.05 while the new one is at 2016.77. While the difference in per data point is very small, as 
is the change in the phase shift parameters, the difference in x^ is important for the calculation of error estimates. 

The procedure used is first to calculate the "natural" phase, (j)o, obtained from the original fit to the experimental 
data. We then calculate a new phase, the "imposed" phase (p = (j)Q — Sait. In the fitting procedure, points from artificial 
data are included which consist of values of the phase to be imposed at a chosen single energy over a restricted range 
oft. The new phase "data" are included in the non-spin-flip amplitude, F{0), only, the spin-flip amplitude being left 
free to have whatever phase the fit prefers. This additional phase "data" is taken to have very small errors to force 
the desired phase. A constant error (0.001 rad) is taken at each of 5 phase data points. In the modified minimization 
process there are two contributions to x^- 



X'^xi + X; (32) 

where Xd is the part of the x^ from the experimental data points and x^ is that coming from the phase points. 

The total (including the phase data) is minimized but the part of x^ of principal interest is that due to the 
experimental data points, Xd- It is found, as anticipated, that the Xd from the true data increases as the imposed 
phase is chosen farther from the original phase, (t>o- 

The result of Xd for fitting the forward phase up to -t=0.32 (GeV/c)^ at a beam momentum of 1.3 GeV/c is shown 
in Fig. [T] An increase in Xd of 4 gives an estimate of the uncertainty of ±0.2 (GeV/c)~^. 

In the upper part of Fig. [2] one sees that the imposed phase is well fit up to 0.4 (GeV/c)^. Attempting to fit to a 
larger range of t did not lead to a better fit to the desired phase but instead the x^ due to the phase portion became 
larger. It appears that the limit has been reached where a modification of the phase linear in t can easily be used. 
The lower part of Fig. [5] shows the limits of ± 0.2 (GeV/c)~^, an uncertainty of around ± 5%. Thus, the phase is 
reasonably well determined. One might believe that thephase of the nuclcon-nucleon amplitude is better determined 
because of the use of OPE in the higher partial waves [131 • 



IV. NUCLEON-NUCLEON AMPLITUDE 



With some confidence that a phase-shift fit to data provides a reasonably reliable determination of the absolute 
phase, we now turn to the nucleon-nucleon phase obtained from the recent fit by Arndt et al. (l^ . The spin average 
is given by 

m4«++|M|++> + < + -|M| + -»^ ^" + ^^- + ^^° . (33) 

The amplitude must also be averaged over the neutron and proton so that the averaged amplitude becomes 

A = ^M{I = 1) + ^M{I = 0) (34) 

In Fig. [3] are shown several properties of the Arndt et al. [l^ amplitudes at 1.8 GeV/c. The value of qr can 
be extracted by averaging the amplitudes or the cross sections. The top panels show the result of the two methods 
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FIG. 3: Properties of the NN amplitudes at 1.8 GeV/c. Panel (a) shows the result of the absolute value of the averaged neutron 
and proton amplitudes. In panel (b) the solid curve shows the result of averaging the cross sections and taking the square root. 
In both (a) and (b) the dash curve shows a plot of the Gaussian approximation with ur — 5.6 (GeV/c)^^ . Panel (c) shows the 
phase relative to the phase at t = (solid curve). Panel (d) shows the absolute phase over an extended range of t. 



compared with the exponential form given in Eq. [T] The panel (c) shows the variation of the phase from the forward 
value. It is seen that a linear approximation in t is reasonable up to —t of about 0.2 (GeV/c)^. Panel (d) shows the 
variation of the total phase. 

Figure [4] shows the variation of the phase for several incident momenta as extracted from Arndt et al. for 
various beam momenta. Also shown are values from Wallace for spin-independent neutron and proton amplitudes. 
Wallace's values were taken from earlier fits by Arndt's group at lower energies. It is seen that for all cases except for 
the very highest beam momenta the phase increases with = —t so that aj < 0. 



V. POTENTIAL MODELS 



In this section we consider potential models for the interaction in an attempt to relate the phase to underlying 
physical parameters. For the Gaussian potential in the eikonal approximation we can make considerable progress 
analytically. We also consider the exact numerical solution for this form of potential as well as that of exponential 
and Woods-Saxon potential forms. 
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FIG. 4: Dependence of the phase on momentum transfer. The curves labeled with a momentum (PLab) are from Arndt et al. 
ia|. The two dashed curves labeled "Wallace" are the phases from neutron and proton spin averaged amplitudes from Ref. 
II. 



A. Gaussian Potentials 

In this section we assume a potential expressed as 

V{r) = Vne-'''/''^ +iV/e-'^'/'''. 
The integral on z in Eq. [3] can be easily done to give 



(35) 



(36) 



where a = i^/Trrj^Vn/v = ia' is purely imaginary and (3 = —^/nriVi/v is real and positive (since Vi must be negative). 
Now the eikonal expression is 



poo 

f{q)^^lk bdbJo{bq)y2 
Jo 1 



^0 1 ^ r 



n=l m— 



^ \ ^n—mpm^—(n—7n--\-rjm)b'^/r'j^ 



(37) 
(38) 
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where rj = rj^/r'j. Following the development in Appendix A the full expansion of the forward amplitude in powers of 
a is given by 



/(O) 



-ikrj 



A(0,l,/3)+^^^Ao(-,l,/?) 



We have introduced a generalization of the incomplete Gamma function with u > 0, 



^ n\{u + n)'^ 



(39) 



(40) 



This function can be computed from its expansion over a large part of its range. Some of its properties, including an 
asymptotic expansion for large final argument, are discussed in Appendix B. Another function, useful when the first 
argument of A is non-zero, is 



,-fe 



— ' n\(u + n)'^' 

n— 



A(u, fc, /?). 



Using Eq. 



]for /(O) we can see that (to lowest order in a) 

,^ + A(i,l,/3) 



p — a 



(41) 



(42) 



A(0,l,/3) 

where A(0, 1, P) < 0. The proportionality of p to a' shows that it is strongly influenced by the real part of the potential. 
The total cross section, in lowest order in a, is given by 

a^2nrj[gi+liiif3)+Ei{p)] (43) 

where gi is Euler's constant. Since, for the range of values of /3 used here, i?i(/3) is small, the dependence of on /3 
is small so that this equation provides a strong constraint on r/. The lowest order correction in a is a^. This is an 
exact representation of the simple eikonal approximation (i.e. without corrections given by Wallace [ll|) for a purely 
absorptive potential. 

Equation IA14I can be expressed as 

■1 2 ^ 



1 



A(0,j + l,/?) + ^ 



^Ao(-,j + l,/3) 



(44) 



We can write the ratio to the forward amplitude (Eq. [55)1 as 

/(g) _ Er=o jr [K0,J + l,/3) + ET=i ^Ao(|, J + l,/3)] (4^)' 



/(O) 



A(0,l,/3) + Efci^Ao(|,l,/3) 



so that, retaining only the first order contribution in a' and q^, we have. 



/(?) 



r,V A(0,2,/3) 



/(O) 4 A(0,l,/3) 

From this expression we can identify 

^1 A(0,2,/3) J 
2 A(0,1,/?) 1 



qr + mj = 



1 + ia 



,Ao(^,l,/3) . ,Ao(l,2,/3) 



A(0,l,/3) 



A(0,2,/3) 



Ao(^,l,/3) Ao(^,2,/3) 
A(0,l,/3) A(0,2,/3) 



(45) 



(46) 



(47) 



Thus, in this linear approximation in a = ia' , we see that is relatively stable under variations in the radii. 
However, the two ratios multiplying a' are very similar for equal radii so that a/ varies rapidly for small variations in 
the radii of the real and imaginary potentials. 

Figure [5] shows several observables obtained from exact solutions for the Gaussian potential of as a function of the 
radius of the real part of the potential for — ±20 MeV. The sign change of p and a/ with sign change of the real 
potential noted in the introduction in the eikonal prescription is approximately reproduced in the exact calculations. 
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FIG. 5: Dependence of several observables on the r.m.s. radius of the real Gaussian potential well for both signs of its strength. 
The solid (dashed) curve shows the result for the positive (negative) potential. Panel (a) shows that the phase parameter has 
a very strong dependence on the radius of the real potential. In panel (b) it is seen that the ratio of real to imaginary part of 
the forward amplitude also has a fairly strong dependence while the Gaussian parameter controlling the fall-off of the absolute 
value and the total and reaction cross sections have little or no dependence on this parameter. 



B. Exponential Potential 



In this section we look at potentials of the form 

V{r) = ly^e-''/'^ + iWie-''/''' . (48) 

We take both Wji and Wj to be functions of the incident momentum in order to fit the data. We have used the phase 
shift analysis of Arndt et al. [l2i | to calculate the spin-isospin average of the amplitude and used these amplitudes as 
a guide to fitting the potentials so we are using a spin-zero on spin-zero calculation to fit spin averaged data. We also 
used the older results of Wallace [l^ as well as the values of p from the Particle Data Group [l^. Figure [H] shows the 
results. 

We see (solid line, direct fit) that the phase parameter a/ is negative at low incident momentum and appears to 
be nearing zero at higher momenta. The parameter p shows a similar behavior with opposite sign as expected from 
the simple (5-function model of Section II. The parameter qr rises slightly at low momenta and the total and reaction 
cross sections are relatively flat although the reaction cross section goes to zero as the momentum is reduced below 1 
GeV/c. 
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FIG. 6: Dependence of several parameters on the radius of the real potential well for the exponential potential. Panel (a) shows 
that the phase parameter has a very strong dependence on the radius of the real potential. The curves are labeled by the real 
radius (in fm). These labels identify the curves for the remaining panels as well. In panel (b) it is seen that the ratio of real to 
imaginary part of the forward amplitude also has a fairly strong dependence while (panels c and d) the parameter controlling 
the fall-off of the absolute value and the total and reaction cross sections have little or no dependence on this parameter at the 
higher energies. The solid triangles denote the values used by Franco and Yin [^. The open circles and squares where taken 
from Wallace p^ . The solid squares were extracted in the present work from the fit of Arndt et al. [l^ and were fitted to get 
the solid curve. The solid dots with errors in panel (b) were taken for the Particle Data Group with several points at low 
laboratory momentum with negative values excluded. 



Also shown in Fig. [5] are curves for changes from the fitted r.m.s. value for the radius of the real potential of 1.21 
fm to 1.33 fm (dot-long dash) and 1.45 fm (dot-short dash) while holding all other parameters fixed.. It is seen that 
there is a significant change in a/ and a corresponding change in p as well. There is also a more moderate (except 
at low momentum) change in aj^. At 2.0 GeV/c aj goes from -2.01 to -3.52 to -5.70 (GeV/c)~^ as the percentage 
change in radius goes from to 10% to 20%. 
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FIG. 7: Dependence of several parameters on the radius of the real potential well for the Woods-Saxon potential. Panel (a) 
again shows that the phase parameter has a very strong dependence on the radius of the real potential. The curves are labeled 
in panel (a) by the rms value of the radius of the real Woods-Saxon potential (in fm). These labels identify the curves for the 
remaining panels as well. The curves and symbols have the same meaning as in Fig. [S] 



C. Woods-Saxon Potential 



In this section we look at potentials of the form 

V{r) = . (49) 

We do not believe that this form provides a realistic representation of the distribution of the strength of the interaction 
for the nucleon-nucleon system but we include it to show that the variation with percent increase in the radius of the 
real potential is about the same as the Gaussian or exponential potential. 

Figure [7] shows results for a fit with this potential. We see that while the basic fit has a larger r.m.s. radius for 
the real potential than for the exponential potential, the change with percentage change in radius is qualitatively very 
similar. At 2.0 GeV/c a/ goes from -1.96 to -3.13 to -4.68 (GeV/c)^^ as the percentage change in radius goes from 
to 10% to 20%. 
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VI. MULTIPLE SCATTERING RESULTS 



Since multiple scattering of the projectile on a nuclear target depends on the phase of the amplitude it may be 
possible to measure the phase (including the phase variation) in this manner. Of course, one must be aware that what 
is obtained is the effective value in a nuclear medium. 



A. Light Nuclei 



Franco and Yin [3,] considered the amplitude as defined in the introduction where it was assumed that only the 
non-spin-flip amplitude was needed for the elastic scattering on light nuclei that they considered and the variation of 
phase was relative to the phase at zero degrees. We follow this same definition. They pointed out that if the ratio 
of real to imaginary part of the amplitude were zero, in the eikonal approximation that they used, the cross section 
would be independent of the sign of a/ . Since this ratio is small, but non-zero, this symmetry is only approximate and 
they found about the same large improvement in the agreement with data [20] for values of -1-10 and -15 (GeV/c)~^ 
at 1.75 GeV/c, per nucleon, -^7.5 and -13 (GcV/c)-^ at 1.25 GeV/c per nucleon and +11.5 and -12 (GeV/c)"^ at 
1.08 GeV/c per nucleon. Clearly the values from the phase shift analyses favor the negative sign although the trend 
with incident momentum seems to be contrary to that seen in Ref. However, except for the 1.75 GeV/A case, the 
values of are very different from those obtained from spin-averaged amplitudes (see Fig. [6] or [7]). A second way 
to obtain the slope parameter is to fit the differential cross section with an exponential in t. In this case the spin-flip 
cross section is included which leads to a more nearly isotropic cross section and hence smaller values of ajj. Thus, 
with this choice some effect of the spin flip is included in the representation. The values for an are believed to be 
equivalent at about 600 MeV [21 1. 

El-Gogary et al. [3] also treated a — a scattering and found (using different NN amplitude parameters) that a phase 
factor linear in the momentum squared of a value of -1-5 (GeV/c)^^ greatly improved the agreement with the data at 
1.75 GeV/c/A. 

Usmani et al. [7] considered the case of a modified helium wave function consisting of the sum of two Gaussian 
pieces. They found that for a scattering on *He only moderate corrections were seen. Since the data extend to a value 
of of 4 (GeV/c)^ this may seem to be contrary to expectations. However, the a'' He scattering is dominated by a 
large number of scatterings in this momentum transfer range so that a typical value of for one of the scatterings 
will be reduced by a factor equal to the number of scatterings, typically between 8 and 16. So the form factor needs 
to be accurate only up to a range of the order of —t = 0.25 — 0.5 (GeV/c)^. 



B. Heavier Nuclei 

Treating heavier nuclei, Lombard and Maillet Q considered, in addition to the cross section, the asymmetry. A, 
and the spin rotation parameter, Q so they were forced to include the spin dependent amplitude. While they used 
for the spin-independent amplitude a basic form with no phase variation, the basic form of the spin-flip amplitude 
included a phase variation consistent with phase shift analyses [see their formula (2)]. The phase variation in their 
basic spin-flip amplitude can be represented by a coefficient of t of about 0.95 (GeV/c)"'^. They then introduced a 
global phase variation which was applied to both amplitudes and estimated the phase parameter which would improve 
the agreement with the cross section and A (Q had not been measured at the time). They found that while A was not 
very sensitive to a variation in this phase, the quantity Q was. However, since there seems to be no reason to believe 
that the two amplitudes would have the same phase variation (or the difference they used), one can only conclude 
from their results that there is a sensitivity to the phase of the amplitudes. 

Lassaut, Lombard and Van de Wiele [5| considered an additional phase variation to be applied to representations 
of both amplitudes determined from amplitude analyses [q^]. They also applied the same phase variation to both 
amplitudes. If one assumes that the non-spin-flip amplitude dominates the elastic cross section and adds the phase 
parameter they obtained (-0.25 fm^=-6.4 (GeV/c)~^) to the one present in their basic amplitudes (-5.6 (GeV/c)~^) 
one arrives at a value of (-12 (GeV/c)~^) for their incident momentum of 1.46 GeV/c, similar to the negative values 
-13 and -15 (GeV/c)~^ obtained by Franco and Yin at nearby incident energies. 

Auger and Lazard considered the effect of multiplying a global phase times a parameterization of both of the 
nucleon amplitudes. Since they considered only asymmetry and spin rotation observables it is difhcult to compare 
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FIG. 8: Comparison of the cross section for aa scattering calculated through various orders for aj = and p = —0.23. The 
data are from Satta et aZ.l20l 



with their results. They considered phase factors with parameters close to those of Franco and Yin, however, since 
these must be added to the effective values already implicit in the phenomenological amplitudes they are not really 
comparable. 

Chaumeaux et aZ.(2^ also considered the addition of a phase variation in the scattering from heavier nuclei. 

One expects that for all but the very lightest nuclei the sensitivity to such a phase will be small since the scattering 
at high energies is mainly determined by the radius and the diffuseness of the surface l^sj. The influence of such a 
phase would mainly be contained in the region of the minima where many different multiple scattering and medium 
effects contribute. 



C. Extracting the Phase Variation from Multiple Scattering 



We have revisited the Glauber calculation of Franco and Yin JS] with a view to determine the effect of the possible 
nuclear modifications on the extraction of a/ from the data [20i | at 1.75 GeV/c/A beam momentum. For helium 
scattering from helium it is possible to do the full multiple scattering calculation in the eikonal approximation. An 
additional advantage is that the basic form factor of helium has no zero in the momentum range of interest whereas for 
heavier nuclei that is not true. Thus, the minima are better understood. They are not, however, simple interferences 
between different scattering orders (as they are approximately for proton scattering on ^He) but arise from a more 
complicated set of interferences. Figure [5] shows the calculation with the amplitudes obtained from the sum up to a 
limited number of orders for a/ = and p — —0.23. One sees, for example, that the first minimum is well defined 
in position only by summing through eighth order scattering and its depth is established by twelth order. One also 
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FIG. 9: Dependence of the cross section for aa scattering on p and aj. The lower panel shows that for p — 0.20 the value of 
ai = —9 (GeV/c)~^ gives an adequate fit to the data of Satta et aZ. [20|. 

sees that the first order scattering is negligible above about 0.6 (GeV/c)^. For p — aj ^ 0, the individual orders of 
scattering are purely imaginary and alternate in sign. The effect of a finite value of p is to multiply the nth order 
amplitude by (1 — ip)". 

One can be concerned that the effects of p and a/ may be confused. The result for the full multiple scattering 
calculation for aa scattering is weW below the data for large while in the case of the heavy nuclei, this is not in 
general true. There the major effect of setting a/ to a non-zero value is seen in the minima. Figure [9] shows results 
for variations in p over the uncertainty observed in the NN data. The major effect for a/ = is to be seen only in the 
minima. For a finite value of a/ of -13 (GeV/c)^^, this is not the case with the cross section away from the minima 
being increased as well. 

In an attempt to see how well the values of a/ are determined from the a ^He data we have calculated a measure 
with regard to the data by Satta et al. [23|. Since no tabulated data were given we took the data from the plots 
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including 20% errors for most points. Figure [TUl shows x^/N for four selected values of p where it is seen that there are 
two minima. If p were zero, the /N curve would be symmetric. For p = +0.20 the value of about -10 (GeV/c)~^ 
is favored for the negative solution. 

Looking at Figs. [S] and [7] it is seen that these values of p and a/j might be possible by a combination of increasing 
the radius of the real part of the interaction and lowering the energy for the evaluation of the NN parameters although, 
even with these assumptions, one is at the very limit. 

Before one can consider such a reconciliation of the multiple scattering and free space determinations of the phase 
variation there are a number of corrections which must be treated. 

a) The spin-flip may give substantial corrections through double spin flip. 

b) The off-shell corrections could give a contribution to the amplitude of similar nature to that of p and ajj. 

c) Short-range correlations in the "^He wave function will modify the form factor. 

d) Three-body forces could affect the scattering. 

e) Standard corrections to Glauber theory [ll| still need to be considered. 
This list is not meant to be exhaustive. 




:p=+0.20, +0.10, -0.10, -0.23 a^=5.6 (GeV/c)-^ : 
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-20 -15 -10 -5 5 10 15 20 

aj (GeV/c)-^ 

FIG. 10: Values of from a comparison with the measured cross section for aa scattering [53] as a function of ai for four 
different values of p. 



VII. CONCLUSIONS 

It has been shown that the phase of the strong scattering amplitude is a sensitive function of the relative size of the 
radii of the real and imaginary parts of a potential describing the interaction. This was done with the consideration 
of three models of the potential. 

We have found that the phase parameter is moderately well determined for the example case of K+p scattering 
from standard amplitude analyses with the controlling principle being the unitary expansion of the amplitude. It is 
expected that the phase in the nucleon-nucleon case would be better determined. 
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Since the effective amplitude in the nucleus can be taken as the free amplitude evaluated at a shifted (normally 
lower) energy one can perhaps understand the larger magnitude of the phase parameter than that seen in free space. 
Studies of this shift in the effective energy of the scattering amplitude have been made in the case of pion-nucleus 
scattering |24j] and the predictions of such an energy shift were verified psj experimentally. Similar corrections have 
been calculated in this energy range for nucleon- nucleus scattering (26| . 

This explanation may not be adequate to give a fit to the data so one is led to consider the possibility of a larger 
radius for the real part of the interaction due to a partial dcconfinement in the nuclear medium. If this is true it may 
be a new way to study "non-classical" modifications in the nuclear medium. 

This work was supported by the National Science Foundation under Contract PHY-0099729. 



APPENDIX A: EXPANSION OF THE NN AMPLITUDE 



Starting from Eq. [38] and using 



we can write 
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This equation is a generalization of Eq. 4.11 in Wallace |ll| for two independent radii. In the forward direction 
(q=0) we can write 
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Transforming the variable of integration from t to z with z = e we have 



/(O) = ^ / dz'- ^ 1. 

Jo ^ 

Since a is thought of as smaller than (3 (in absolute magnitude) we expand the exponential in a to find 

no) f dz. 
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Since 



1 I -f3z _ il oo (_r,^n 
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^ 71—1 



(A9) 



where gi is Euler's constant (= 0.57721 . . .) and Ei{p) is the exponential integral, we have a closed form for the 
leading order which leads directly to Eq. 03] in the main text. 
The full amplitude reads 



dt 



and upon using the same change of variable 
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Expanding the a exponential we have 
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and since 
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we may finally write 
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APPENDIX B: PROPERTIES OF THE FUNCTION A 

The power series definition of A(u, k, x) is convergent for all values of u and x but for large values of x (greater 
than about 40) it is numerically difficult to calculate in that manner. We note that the relation 



x^^%^=A(0,fc-l,:.) 
dx 

follows from the series definition. From this equation and from the known relation for fc = and fc = 1 

A(0, 0, x) = e-"" - 1; A(0, 1, a;) = - ln(a;) - gi - Ei (x) 

we can iterate to find the asymptotic behavior for large x. 
The asymptotic behavior of A(0, k, x) is given by 



(Bl) 



(B2) 



A(0,fc,a:)^-^ 



gk-t[\n{x)Y 



(B3) 



for fc > 0. Note that A(0, fc, x) < and A(0, fc, a;) > A(0, k + l,x) for x > and aU fc > 0. 
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1.000000 


1 


0.577215664 


2 


0.9890560 


3 


0.9074791 


4 


0.9817279 


5 


0.9819955 



TABLE I: Values of the coefficients gj determined numerically (except for the first two). 



The values of X{0,k,x) can be calculated to 5 significant figures with the constants given in Table |T] for x greater 
than 10. The values for x less than 40 can be easily calculated by the series. 
For the variation in the first parameter, u, the relation 

Xiu + l,k,x) = ^-^-—^ (B4) 

allows the calculation the function for all values of u from those between and 1. For small u we can expand the sum 
to get 

A(u, k, x) = A(0, k, x) - fcuA(0, k+l,x) + ... (B5) 
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